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Continental runoff is a major source of freshwater, nutrients and terrigenous material to the Arctic Ocean. 
As such, it influences water column stratification, light attenuation, surface heating, gas exchange, biological 
productivity and carbon sequestration. Increasing river discharge and thawing permafrost suggest that the 
impacts of continental runoff on these processes are changing. Here, a new optical proxy was developed and 
implemented with remote sensing to determine the first pan- Arctic distribution of terrigenous dissolved 
organic matter (tDOM) and continental runoff in the surface Arctic Ocean. Retrospective analyses revealed 
connections between the routing of North American runoff and the recent freshening of the Canada Basin, 
and indicated a correspondence between climate-driven changes in river discharge and tDOM inventories in 
the Kara Sea. By facilitating the real-time, synoptic monitoring of tDOM and freshwater runoff in surface 
polar waters, this novel approach will help understand the manifestations of climate change in this remote 
region. 
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Continental runoff influences surface ocean processes throughout the Arctic. Rivers discharge about 
3,300 km 3 yr" 1 of buoyant freshwater to the Arctic Ocean 1 where it sustains a halocline that isolates 
the upper 200-m layer from a warmer, saltier Atlantic layer 2 . This upper layer represents only 0.1% of the 
global ocean volume but receives 11% of global riverine discharge from a circumpolar drainage basin that 
encompasses over 10% of the world's land surface area. Riverine input is evident throughout the Arctic from 
low salinities and high concentrations of terrigenous dissolved organic matter (tDOM) in polar surface waters 3 ' 4 . 
Arctic rivers discharge onto the world's largest continental shelves that surround and cover more than half of the 
Arctic Ocean area 5 . These river-influenced ocean margins are dominant features of the Arctic and important 
regions of sea-ice production, biological productivity, carbon sequestration and processes that drive biogeochem- 
ical cycles on a basin-wide scale. 

The impact of rivers on the freshwater content and biogeochemistry of the Arctic Ocean is strongly influenced 
by the residence time of riverine inputs and their pathways to the North Atlantic Ocean. The average residence 
time of Arctic shelf water ranges from months to years depending on shelf region 6 ' 7 . Variations in atmospheric 
forcing can affect residence time by dramatically altering continental runoff routing and distribution in the Arctic 
Ocean 7,8 . Climate-driven changes in runoff distribution have been linked to large inter-annual and decadal 
variations in the freshwater content of the Arctic Ocean, and these fluctuations have important ramifications 
for North Atlantic meridional overturning circulation 9 . Furthermore, changes in runoff routing through the 
Arctic Ocean can regulate the extent to which Arctic tDOM is incorporated into North Atlantic Deep Water 
(NADW) and distributed in the global ocean 4 . Despite their global significance, continental runoff routing and 
distribution in the Arctic remain poorly understood and difficult to predict in this remote and under-sampled 
region of the global ocean. Here, a new optical proxy for tDOM implemented using remote sensing of ocean color 
facilitates the synoptic and retrospective monitoring of continental runoff in the Arctic Ocean, and provides a 
novel approach to assess climate-driven changes in the Arctic. 
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Results 

Dissolved lignin is well established as a biomarker of tDOM in 
oceanic waters 1011 and has been successfully applied as a tracer of 
riverine inputs in the Arctic Ocean 3 ' 4 . However, its use relies on the 
collection of seawater samples and land-based laboratory analyses 
that prevent synoptic applications. Lignin is also an important chro- 
mophore in tDOM, a property that facilitates detection using optical 
properties. Here, we demonstrate that the spectral slope coefficient of 
chromophoric dissolved organic matter (CDOM) between 275 and 
295 nm (S275-295) is a practical and reliable optical proxy for dis- 
solved lignin and tDOM across various river-influenced ocean mar- 
gins of the Arctic Ocean 1213 (Fig. la). An increase in S275— 295 is 
indicative of a decrease in the tDOM content of surface waters. 
Under the assumption that S275-295 behaves conservatively, an 
increase in $275-295 corresponds to a diminishing influence of con- 
tinental runoff. Biological degradation has little effect on S275-295 13 ' 14 
and photodegradation of tDOM in polar surface waters is not extens- 
ive 415 , thereby warranting the assumption and the use of S275-295 as a 
tracer of tDOM and freshwater runoff in river-influenced ocean 
margins of the Arctic. 

The practicality of the S275-295 proxy resides in its amenability to 
remote sensing. The development of an empirical algorithm reveals 
that S 2 75_295 can be estimated from ocean color with excellent and 
consistent accuracy (within ± 4%) over the broad range of water 
types and S275-295 values observed in the Arctic Ocean (Fig. lb). The 
strong connection between S275-295 an d ocean color in river-influ- 
enced ocean margins exists because the optical property $275-295 is an 
excellent tracer of tDOM 13 and the gradients of ocean color in river- 
influenced ocean margins are primarily driven by the constituents of 
continental runoff. The S275-295 algorithm thus facilitates the syn- 
optic and real-time monitoring of tDOM and freshwater runoff in 
the Arctic using ocean-color remote sensing. 

Implementation of the S275-295 algorithm using MODIS Aqua 
ocean color provides the first pan-Arctic view of average tDOM 



distributions in river- influenced ocean margins during August 
2002-2009 (Fig. lc). The influence of the three largest Arctic rivers 
on tDOM distributions is evident along the entire Siberian margin. 
One of the most striking features is the overwhelming amount of 
tDOM in the Laptev and East Siberian Seas, indicating a large portion 
of Eurasian river tDOM and runoff is routed alongshore in an east- 
ward direction. This observation is supported by measurements of 
high in situ concentrations of dissolved organic carbon of terrigenous 
origin in the East Siberian Sea in August 2008 16 , and is consistent 
with shelf water residence times of several years in this region 6 ' 7 . Low 
S 2 75_295 values along the Alaskan coast in the Chukchi Sea depict the 
transport of Yukon River tDOM to the Arctic Ocean in the Alaska 
Coastal Current and indicate notable tDOM contributions by smal- 
ler Alaskan rivers (e.g., Noatak and Kobuk Rivers). Most of the runoff 
and tDOM from the Mackenzie River enters the Beaufort Sea and 
appears to be entrained in the gyre circulation of the southern 
Canada Basin. A relatively small fraction of the discharge from the 
Mackenzie River was routed to the Labrador Sea via the Canadian 
Arctic Archipelago (CAA). 

Assessment of inter- annual variability in tDOM distribution 
revealed a change in the routing of Mackenzie River discharge over 
the last decade (Fig. 2; Supplementary Figs. S3 and S4). The 
Mackenzie River outflow progressed from an alongshore, eastward 
path to the CAA in 2002 (Fig. 2a) to a cross-shelf, northwestward 
path to the Canada Basin since 2006 (Fig. 2b). This routing switch has 
important implications for the fate of North American runoff and 
tDOM. Routing through the CAA favors the export of North 
American tDOM to the Labrador Sea 17 . Although the magnitude of 
this export remains unknown, deep ocean convection in the 
Labrador Sea can sequester tDOM in North Atlantic Deep Water 
for centuries 4 . In contrast, routing to the Canada Basin contributes to 
longer-term storage of runoff and processing of tDOM within the 
Beaufort Gyre and Arctic Ocean 4 ' 15 . This decadal shift to greater 
export of Mackenzie River runoff to the Canada Basin during 
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Figure 1 | The CDOM spectral slope coefficient, S 2 75-295> as a tracer of riverine inputs in the Arctic Ocean, (a) The relationship between 
S275-295 and dissolved lignin carbon yield (TDLP 9 -C) across various Arctic river-influenced ocean margins, (b) A performance evaluation of the S275-295 
algorithm. On average, S275-295 can ^ e estimated from ocean color within 4% of S275-295 va l ues measured in situ, (c) Implementation of the 
algorithm using MODIS Aqua ocean color providing a pan-Arctic view of an August climatology (2002-2009) of £275-295- Increasing S275-295 values are 
indicative of a decreasing fraction of tDOM. The four largest Arctic rivers are labeled and ranked in order of decreasing discharge: Yenisei (1), Lena (2), 
Ob (3), and Mackenzie (4). River-influenced margins of the Arctic are labeled: Gulf of Ob (GO), Kara Sea (KS), Laptev Sea (LS), East Siberian Sea 
(ESS), Chukchi Sea (CS), Beaufort Sea (BS) and Amundsen Gulf (AG). The contour lines represent the 2000-m isobath. 



SCIENTIFIC REPORTS | 3 : 1053 | DOI: 1 0.1 038/srep01 053 



2 




0.015 0.020 0.025 0.030 0.035 



Figure 2 | Inter-annual changes in the routing of Mackenzie River runoff, (a) The Mackenzie River outflow in 2002, showing eastward transport of 
North American runoff to the Canadian Arctic Archipelago, (b) The Mackenzie River outflow in 2008, showing northwestward transport of North 
American runoff to the Canada Basin. A progressive switch from eastward to northwestward routing occurred between 2002 and 2011 (Supplementary 
Figs. S3 and S4) and coincides with the rapid freshening of the Canada Basin. The $275-295 algorithm was implemented using 4-km resolution, yearly- 
binned, MODIS Aqua ocean color. The Mackenzie River is labeled (4). The contour line represents the 2000-m isobath and outlines the Canada Basin. 



2006-2011 could contribute to less North American tDOM being 
sequestered in the deep ocean. 

The routing of Mackenzie River outflow since 2006 coincides with 
observations of a rapid accumulation of freshwater in the Canada 
Basin. Freshening began in the 1990s and accelerated in the late 
2000s 18 " 21 . The realization that a large accumulation of freshwater 
in the Canada Basin could impact global ocean circulation stimulated 
research to identify the freshwater sources and climatic forcing 
responsible for observed changes in salinity. Pacific water, ice-melt, 
precipitation and river runoff are distinct sources of freshwater to the 
Canada Basin. The relative contributions of these freshwater sources 
to the freshening are typically estimated using tracers (e.g. (5 18 0, 
alkalinity, salinity) or are indirectly calculated using complex 
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models 9,20 ' 22 " 25 . The lignin-based $275-295 P rox Y presented here sub- 
stantiates the riverine origin of freshwater and provides a direct 
means to distinguish continental runoff from other freshwater 
sources using remote sensing. 

The acceleration of the Beaufort Gyre due to the strengthening of 
the Beaufort High 26 and routing of Eurasian river runoff to the 
Canada Basin driven by the Arctic Oscillation 9 are complementary 
and currently leading explanations for the freshening of the Basin 27 . 
The remotely sensed runoff distributions presented in this study 
(Fig. 2 and Supplementary Figs. S3 and S4) further suggest the recur- 
rent influx of Mackenzie River runoff to the Canada Basin in the late 
2000s was also a significant source of freshwater to the Basin and 
contributed to the freshening. The temporal sequence of the 
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Figure 3 | Temporal sequence of the Mackenzie River freshet release into the Canada Basin during June and July 2008. The $275-295 algorithm was 
implemented using 4-km resolution, 8-day-binned, MODIS Aqua ocean color acquired between June 1 and August 3. The discharge value (km 3 ) 
displayed on each panel indicates the cumulative Mackenzie River discharge released between May 1 (start of the freshet) and the end of each 8-day period. 
The contour line represents the 2000-m isobath and delineates the Canada Basin. The Mackenzie River is labeled (4). 
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Figure 4 | Climatic forcing of major Eurasian riverine inputs, (a) The North Atlantic Oscillation (NAO) and Arctic Oscillation (AO) wintertime indices 
(DJFM: Dec-Jan-Feb-Mar) in relation to the combined discharge anomaly of the Yenisei, Lena and Ob rivers between 2002 and 2009. High NAO and AO 
indices are typically associated with higher precipitation over Northern Eurasia. River discharge is calculated as the cumulative discharge from May 1 to 
August 30, and the anomaly is relative to the 2002-2009 average, (b) Average S 275 _ 2 95 values for surface waters of the Kara Sea in relation to the combined 
discharge of the Yenisei and Ob rivers between 2002 and 2009. The S 275 _ 2 95 values for the Kara Sea were calculated as the average value in the region 
displayed in panel (c). River discharge is calculated as the cumulative discharge from May 1 to August 30. (c) Distributions of riverine inputs during high 
(2002) and low (2005) discharge years. The years correspond to low and high average S 275 _ 2 95 values observed in the Kara Sea between 2002 and 2009. The 
^275-295 algorithm was implemented using MODIS Aqua ocean color data averaged from May to October. The Yenisei and Ob rivers are labeled ( 1 ) and 
(3), respectively. 



Mackenzie River discharge in 2008 provides unambiguous evidence 
that the spring freshet was directly and rapidly injected into the 
Canada Basin during June and July (Fig. 3). In situ measurements 
of salinity, c5 ls O, and alkalinity from an independent study 20 corrob- 
orate that Mackenzie River runoff was a minor source of freshwater 
to the Canada Basin in the early 2000s but increased in the southern 
part of the Basin in 2007. 

The contribution of North American runoff to the Canada Basin 
freshwater budget appears to be shifting in response to climatic for- 
cing. A freshwater budget based on data from 2003-2004 estimated 
about 40% of North American runoff enters the Canada Basin 22 . The 
mean freshwater residence time in the Basin is about 10 years 22 , so 
this estimate is representative of the decade leading up to 2004. 
Although significant freshwater contributions from the Mackenzie 
River to the Basin were observed during that decade 28 , this estimate 
does not integrate the increasing contributions from North Ameri- 
can runoff since 2006. The detailed progression of the 2008 
Mackenzie River outflow indicated up to 155 km 3 of Mackenzie 
freshwater was released into the Canada Basin in June and July alone 
(Fig. 3), indicating more than half of the discharge from the Mac- 
kenzie has been entering the Canada Basin since 2006. Increasing 
contributions of North American runoff to the Canada Basin are 



consistent with enhanced Ekman pumping resulting from the 
wind-driven spin-up of the Beaufort Gyre 26 and therefore appear 
to result from recent trends in climate variability. 

Discharge from major Arctic rivers has increased over the past 
decades in response to climate variability 29 " 31 . Climatic forcing of 
river discharge is evident between 2002 and 2009 from the strong 
correspondence between the North Atlantic Oscillation (NAO) or 
Arctic Oscillation (AO) wintertime indices and the river discharge 
anomaly of Eurasian rivers (Fig. 4a). Manifestations of inter- 
annual discharge fluctuations in the Arctic Ocean are best illu- 
strated in the Kara Sea, which receives inputs from the Ob and 
Yenisei rivers. The Kara Sea is a semi-enclosed margin that is 
largely ice-free during summer, which facilitates inter-annual 
comparisons of remotely sensed variables. Inter-annual compar- 
isons revealed consistency between discharge from the Ob and 
Yenisei rivers and the average S275-295 value for the Kara-Sea 
sector (Fig. 4b). The correspondence between discharge and 
tDOM in the Kara Sea is particularly striking when comparing 
2002 and 2005, which are representative of high and low discharge 
years between 2002 and 2009 (Fig. 4c). This demonstrates climate- 
driven changes in the tDOM inventory of the surface Arctic 
Ocean are traceable using this approach. 
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Discussion 

The S275-295 P rox y an d algorithm presented in this study facilitate the 
real-time and synoptic monitoring of tDOM and freshwater runoff 
in surface polar waters. The foundation of the $275-295 P rox y on the 
lignin biomarker is fundamental to the approach. The "terrigenous" 
attribute of lignin substantiates the riverine origin of dissolved 
organic matter and freshwater. Its "chromophoric" property, on 
the other hand, facilitates implementation with remote sensing. 
Although the Arctic remains one of the most challenging environ- 
ments for remote sensing, it is also a region where remote sensing has 
much to offer given the logistical restrictions and high costs of field 
surveys, particularly on the Siberian shelves. The applicability of the 
approach will improve as the extent and duration of sea-ice cover 
decrease, allowing for greater spatial and temporal coverage. 

New capabilities to monitor surface tDOM distributions will prove 
valuable for understanding future climate-driven changes in the bio- 
geochemistry of the Arctic. Permafrost thawing and precipitation in 
drainage basins of Arctic rivers are projected to intensify with escal- 
ating atmospheric temperatures and thereby enhance the mobiliza- 
tion of soil organic matter to the Arctic Ocean 32 36 . The fate of this 
material and its effects on biogeochemical cycles will depend on its 
routing and residence time in polar surface waters 4 . Furthermore, 
processing of organic matter in surface waters is influenced by factors 
that are currently being altered by climate change (e.g., ice cover, 
water temperature). Multiyear records of remotely sensed tDOM 
distributions provide direct evidence of how the routing, inventory, 
storage and residence time of tDOM in surface polar waters change 
in response to climatic forcing. Such applications will thus provide 
crucial information for understanding the fate of immobilized soil 
organic matter and its impacts on biogeochemical cycles in this 
rapidly changing Arctic environment. 

Finally, remote sensing of continental runoff provides useful 
insights about the sources of freshwater in the Arctic Ocean. Ice melt, 
precipitation and runoff are all increasing under the current climatic 
trends 31 ' 37 , and are altering the freshwater budget of the Arctic Ocean. 
Remotely sensed runoff distributions provide direct evidence of con- 
tinental runoff contributions to specific regions of the surface Arctic 
Ocean, and will thereby help understand the mechanisms respons- 
ible for future changes in the Arctic freshwater budget. It is important 
to remember, however, that the uncoupling of river water from its 
tDOM component during the winter freeze-thaw cycle 38 hampers the 
use of this approach to assess year-to-year changes in freshwater 
runoff storage. Future applications will also likely include studying 
the influence of continental runoff on biological processes, such as 
primary production, in polar surface waters. 

Methods 

Surface water samples and data (n = 106) used to establish the relationship 
between S275-295 and TDLP 9 -C (Fig. la) were acquired during several field 
expeditions to the Arctic Ocean between 2005 and 2010, including the Amundsen 
Gulf, Gulf of Ob, and the Beaufort, Chukchi, Laptev and East Siberian Seas 
(Supplementary Fig. SI). Samples for DOC and CDOM were filtered onboard and 
analyzed in the home laboratory using high temperature combustion 39 and a dual- 
beam UV-visible spectrophotometer 1213 , respectively. Samples for dissolved lignin 
(1-10 L) were extracted onboard 40 and analyzed using the cupric oxide oxidation 
method and gas chromatography with mass spectrometry 41 . DOC-normalized 
yields of nine lignin phenols (TDLP 9 -C) were calculated 13 , and the non-linear 
regression of TDLP 9 -C on S275-295 utilized the exponential function shown in 
equation (1), 

TDLP 9 -C= exp(a + jSS 27 5- 295)+ exp (yS 275 - 295) + <Sexp(S 275 - 295) (1) 

where a = 2.7428, = -245.5899, y = 1.6669, and d = -0.9932 are the derived 
regression coefficients. 

The in situ measurements of ocean color (n = 236) and concurrent CDOM 
samples used in the S 2 7 5 _295 algorithm development (Fig. lb) were collected during 
the Malina and ICESCAPE field expeditions (Supplementary Fig. S2). Multispectral 
remote- sensing reflectances, R rs {^)> were derived for 5 wavelengths (/I = 443, 488, 
531, 555 and 667 nm) from profiles made with a Biospherical Compact- Optical 
Profiling System (C-OPS) 42 , and the corresponding CDOM samples were collected 
and analyzed according to NASA standard protocols 43 . "Measured" S275-295 were 



calculated from CDOM absorption coefficient spectra, and regressed on R rs (/l) using 
the multiple linear regression described in equation (2), 

ln[S 275-295] = tf + /?ln[-Rrs(443)] + y ln[.Rrs(488)] + 

( 2 ) 

<51n[2to(531)] + eln[£rs(555)] +£ln[£rs(667)] 

where a = -3.4567, (3 = 0.4299, y = 0.0924, S = -1.2649, e = 0.8885, and 
£ = —0.1025 are the derived regression coefficients. This data set covers abroad range 
of water types found in river- influenced ocean margins, including high-CDOM and 
sediment-laden coastal environments, productive waters at intermediate salinities 
and oligotrophic waters at higher salinities. The inclusion of data from the Chukchi 
Sea in summer (ICESCAPE) further demonstrates that the algorithm performs very 
well in areas where tDOM is relatively low and biogenic particles are abundant and 
variable. Similar accuracy of S 2 7 5 _295 retrievals was obtained using data collected 
seasonally in the northern Gulf of Mexico in various water types ranging from 
riverine to oligotrophic, thereby indicating that the approach used in the algorithm is 
robust. 

MODIS Aqua data were accessed from the NASA ocean color website (http:// 
oceancolor.gsfc.nasa.gov/). Monthly NAO and AO indices were acquired from the 
NOAA Earth System Research Laboratory website (http://www.esrl.noaa.gov/psd/ 
data/climateindices). The Mackenzie River discharge data were obtained from the 
Environment Canada Water Office website (http://www.wateroffice.ec.gc.ca). 
Eurasian river discharge was obtained from the ArcticRIMS website (http:// 
rims . unh . edu/ data . shtml) . 
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